Time dependent evolution of RF-generated non-thermal particle distributions in fusion plasmas

We describe fully self-consistent time-dependent simulations of radio frequency (RF) generated ion distributions in the ion cyclotron range of frequencies and RF-generated electron distributions in the lower hybrid range of frequencies using combined Fokker–Planck and full wave electromagnetic field solvers. In each regime, the non-thermal particle distributions have been used in synthetic diagnostic codes to compare with diagnostic measurements from experiment, thus providing validation of the simulation capability. The computational intensive simulations require multiple full wave code runs that iterate with a Fokker–Planck code. We will discuss advanced algorithms that have been implemented to accelerate both the massively parallel full wave simulations as well as the iteration with the distribution code. A vector extrapolation method (Sidi A 2008 Comput. Math. Appl. 56) that permits Jacobian-free acceleration of the traditional fixed point iteration technique is used to reduce the number of iterations needed between the distribution and wave codes to converge to self-consistency. The computational burden of the parallel full wave codes has been reduced by using a more efficient two level parallel decomposition that improves the strong scaling of the codes and reduces the communication overhead.